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Abstract 

We present a method to compute the probability distribution function of the (true) Hubble constant 
and the age of the universe, given the estimate of the Hubble constant in our nearby galaxy samples. 
Our method takes into account both the observational errors and the cosmic variance, and enables to 
quantitatively compute the constraints on the cosmological models. Based on the present local observation 
Ho = 80 ± 17 km/s/Mpc, the probability of Ho < 50 km/s/Mpc is about 6% for the Einstein - de Sitter 
universe (fio = 1) and 4% for an open (fio = 0.2) universe. These probabilities are reduced to 0.8% 
and 0.03%, respectively, if the accuracy of the observational uncertainty is improved within 10%. Similar 
probabilistic constraints on to are also discussed. 

Subject headings: cosmology: theory - large-scale structure of the universe - methods: statistical 
1. Introduction 

The Hubble constant i?o(= lOO/i km/s/Mpc) is one of the fundamental parameters in cosmology which fixes 
the scale of the various important quantities in astrophysics. The recent observations of the Cepheid variables 
in the Virgo cluster by Hubble Space Telescope suggested a somewhat large value, h = 0.8 ± 0.17 (Freedman 
et al. 1994). On the other hand, the Sunyaev-Zel'dovich effect and the time delay due to the gravitational lens 
effect predict smaller values around h ~ 0.5 (e.g., Birkinshaw & Hughes 1994; Yamashita 1994; Gorenstein, 
Falco, & Shapiro 1988; Rhee 1991). Although the latter methods contain considerable uncertainties, these 
estimates may be close to the true (global) value because these are the high redshift observations. In general, 
lower h is desirable in order to reconcile the age of the universe t ~ H^ 1 with that of the oldest stars in the 
globular cluster (e.g., Chaboyer et al. 1992). 

It is possible that the Hubble constant determined from the local observations systematically differs from its 
global value due to the density inhomogeneity on the scale of the observational samples; the expansion rate of 
small patches of the universe varies from place to place depending on the degree of such inhomogeneities (Suto, 
Suginohara, & Inagaki 1995). In addition, the uncertainties intrinsic to the local observations complicate the 
estimates of the true value of Ho- To what extent can we estimate Ho and to on the basis of the value derived 
from the local observations ? To quantitatively discuss the cosmological implications from the local Hubble 
constant estimate, it is important to derive the probability distribution function (PDF) of the global Hubble 
constant. 

Previously, Turner, Cen, & Ostriker (1992, hereafter TCO) considered the difference between the local and 
global Hubble constant directly from the N-body simulation data, and calculated the PDF numerically. Suto, 
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Suginohara, & Inagaki (1995) derived an analytical expression for the lower limit on Ho and the upper limit 
on to with given value of the local Hubble constant. 

In this Letter, we present a method to compute the PDF of the global Hubble constant using the non-linear 
spherical model. Our semi-analytical expression applied to the cold dark matter (CDM) model reproduces the 
numerical result by TCO very well. We discuss the constraints on the cosmological parameters quantitatively 
in the CDM model as an example. 



2. Matching the Local and Global Universes in a Spherical Approximation 

Let us consider the matching of the homogeneous global universe with the local spherical region of radius r and 
mass M , first in the case of a vanishing cosmological constant A. For open universes, the expansion equations 
are given parametrically as (e.g., eqs.[19.11] and [97.21] in Peebles 1980) 

r=|f (coshfl-l), * = J|^_( s inh0-0), (1) 

for the local spherical region, and 

GM . , 1N GM . . , 

a = — — (coshry - l), t = /2 (smhr) - rj), (2) 

for the global universe, where E and —K are some positive constants. For the later discussion, we define the 
following conventional variables with the subscript L denoting those in the local region: 

r 8wGp L _ 3M 

tlL = ~, "L = » u2 > PL = -, v \ 6 ) 

r 6Hi iirr- 5 



h = -, si = — f, P=~^, 4 



p ' H 
With equations (1) and (2), a straightforward calculation yields 



r PL- P c _ H L -H 

6= , b H = — . (5) 



/ sinh 9 — 9\ 2 ( cosh n — 1 \ 3 
6 =\ ^^ „ „„„u a 1 ( 6 ) 



\ sinh rj — rj J \ cosh 9—1 
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sinh 9 ( sinh 9 — 9 \ / cosh rj — 1 \ 
H sinh rj \ sinh tj — rj J \ cosh 9 — 1 J ' 

rj = Arccosh ^— — 1^ . (8) 

Once SI of the global universe is specified, rj is determined by equation (8), and hence bjj is related to 6 through 
the parameter 9 (eqs.[6] and [7]). Although the above results apply for open (SI < 1 and Sl^ < 1) universes, 
the corresponding results for closed (SI > 1 or Sl^ > 1) models are reproduced simply by replacing 9 — ► i9 or 
rj — ► irj. The relation of 6h versus 6 is plotted in Figure 1 (solid lines; cf, Figure 1 of Davis et al. 1980). 

If A = A/(3if 2 ) is not zero, the above procedure cannot be performed analytically. Assuming that A is 
the same everywhere, we repeated the calculation numerically for spatially flat (SI + A = 1) universes. The 
resulting 6-6 jj relation is plotted also in Figure 1 (dotted lines; cf, Figure 2 of Peebles 1984). Figure 1 clearly 
indicates that the 6-6h relation is insensitive to A, especially in the linear (\6\ <C 1) regime. This is simply 
because A acts as a homogeneous density field which hardly affects the dynamics of the perturbed region 
(Lahav et al. 1991). 
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Figure 1: The relation between the density contrast 8 and the "Hubble constant contrast" 6h for £1 = 0.2, 0.5, 
1 and 2. Solid and dotted curves indicate the relation for A = and SI + A = 1 models, respectively. 



For later use, let us define the following quantity / as 



6 = 



3(cosh r) — l)(r) cosh rj — 3 sinh rj + 2rj) 
sinh rj (4(cosh rj — 1) + sinh ry(sinh rj — 3?y)) ' 



(9) 
(10) 



where rj is related to SI by equation (8). The second line follows from equations (6) and (7). This quantity 
is identical to rfln_D/c?lna where D is the growing solution in linear perturbation theory (e.g., eq.[14.8] in 
Peebles 1980). Again /(SI, A) is very insensitive to A; Lahav et al. (1991), for example, approximate / ~ 
SI 06 + ygA(l + ^Sl) which is very weakly dependent on A. Therefore, we use equations (6), (7) and (10) even 
when A is not zero. 



3. Probability Distribution Function (PDF) of the Global Hubble Constant 



In the above spherical model, 6 and 6h have one-to-one correspondence. Therefore the conditional PDF of 
the global Hubble constant Ph\s{h\6) for given 6, is related to that of the local Hubble constant Ph L {hi,) as 
Ph\6 (h\8)dh = Ph L {hL)dhL- Then the PDF of the global Hubble constant Ph(h) is obtained by the following 
convolution: 



Ph(h) = I P h \ s (h\6)P s (6)d6 

(l + 8 H )P hL (h(l + 8 H ))P s (8)d8. (11) 
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Table 1: Model parameters 



Model 


fio 


A 


hobs 


A/i j s 


^(/i-'Mpc) 


b 


El 


1 





0.8 


0.08 


12 


1 


E2 


1 





0.8 


0.17 


12 


1 





0.2 





0.8 


0.08 


12 


1 


LI 


0.2 


0.8 


0.8 


0.08 


12 


1 


L2 


0.2 


0.8 


0.8 


0.17 


12 


1 


L3 


0.2 


0.8 


0.8 


0.17 


70 


1 


TCO 


1 





0.8 


0.08 


30 


1.3 



Since Ph L {hL) is unknown a priori, we adopt a natural and reasonable choice of Gaussian: 



Ph L (h L ) 



1 

2tt~A/i 



■ exp 



obs 



(12) 



At present, the uncertainty of the observation Ah b s is dominated by systematic errors (see Table 1 of Freedman 
et al. 1994). As the number of the sample galaxies increases, however, statistical errors will dominate and the 
assumption of Gaussian will be more realistic. According to the widely accepted view, the density fluctuation 
field is assumed to be random-Gaussian with the rms value a: 



1 



27T(7 



exp 



"2<7 2 



(13) 



Equation (11) together with equations (6), (7), (12) and (13) yields the analytical expression for the PDF of 
the global Hubble constant. TCO computed this PDF numerically by the similar convolution assuming that 
Ps H \h(^H\h) is the same as P Sh \ hL (6 H \h L ). 

The value of a is computed from the power spectrum P(k) of the density fluctuation: 



1 

2^2 



P(k)W 2 (kR)k 2 dk, W(x) 



-(sin x — x cos x), 



(14) 



where we assumed the top-hat window function W(x) with R being the smoothing length. In what follows, 
we consider the CDM models for specific examples, and adopt the fit by Bardeen et al. (1986) for P(k). The 
value of a is normalized to be & _1 at R = 8/i -1 Mpc (b is the biasing parameter). Note, however, that the 
above result depends on cosmological scenarios only through a (via P(k)). Thus our conclusions below apply 
for any scenarios as long as a takes the same value. 

Our model parameters are summarized in Table 1. We adopt a set of parameters used by TCO for 
comparison. We consider the smoothing length R = 12/i _1 Mpc and 70/i _1 Mpc corresponding to the recession 
velocities of the Virgo (= 1180 ± 22 km/s, Sandage & Tamman 1990) and of the Coma (= 6931 ± 45 km/s, 
Aaronson et al. 1986), respectively. We have numerically integrated the PDF (11) with equations (6), (7), 
(12), (13) and (14), and the results are plotted in Figure 2a. In practice, the integration is performed over 
— 1 < 8 < oo (viz., oo — ► and —> 2iri for 9). For comparison, the Gaussian distribution of P\ lL with 
(h bs , Ah bs) = (0.8,0.08) is shown in dotted curve; the Gaussian curve with (0.8,0.17) is almost identical to 
the model L3, because a is so small that the probability distribution is dominated by the observational error 
(see eq.[17]). The peaks of Ph(h) are slightly shifted to the left of h b s , especially in models El and E2. This 
shift is originated from the convolution process and can be understood by differentiating equation (11) with 
respect to h. 

Since we do not take into account the negligibly small A-dependence in the b-bu relation (Figure 1), the 
results for models O and LI are identical. To see the parameter dependence of the shape of the PDF, we 
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approximate equation (11) using linear theory, i.e., 6h = — /6/3, where / is defined in equation (10): 



Ph(h) 



1 Ah 2 obs +f 2 a 2 h obs h/9 
^(Ah 2 obs +pa 2 h 2 / 9) 3 / 2 



exp 



(h 



l-ob. 



2Ah 2 obs +pa 2 h 2 /9 



(15) 



where integration was carried out formally from — oo to oo with respect to 6. This linear approximation is 
justified only when a is much less than unity and when the integration in equation (11) is contributed mainly 
from the linear regime. Nevertheless equation (15) provides a good insight into the general feature of P^h). 
Equation (15) indicates that Ph(h) is not symmetric around h obs and the tail to the right falls off less rapidly 
than to the left. Hence the expectation value (h): 



(h) 



hP h (h)dh 



Ph(h)dh 



(16) 



becomes slightly greater than h obs even though the peak position is smaller than h obs . It should be noted 
that Ph(h) in the present case is not normalized in an exact sense; the Gaussian distribution (13) assumes 
that 6 ranges from — oo to oo despite that 6 > —1 in reality. In addition, we consider the expanding (h > 0) 
models only. For these two reasons, we define (h) by equation (16). The standard deviation is very roughly 
approximated as 



Ah ee Vdh - W) 2 ) ~ ^K bs + fW ols /9 



Ah 



obs 



)1.2, 

'0 



obs 



l8Ah ob 



(17) 



implying that the width Ah obs of the original PDF is broadened to Ah by the convolution. The effect of 
the convolution Ah — Ah obs increases monotonically with fio as seen in Figure 2a, and decreases with R and 
b simply because a is generally a decreasing function of R and b. In Figure 2b is plotted the cumulative 
probability Ph<{h) that the Hubble constant is less than h, defined by 



Ph<(h) 



Ph(h)dh 



P h (h)dh. 



(18) 



Similarly one can derive the PDF of the age of the universe to or r 
Hq, through P T (r)dT = Ph(h)dh. Combining with 



to /Gyr given the local estimate of 



T 

H~o' 



T(tt ,\o) 



fin 



+ (l-fi - A ) + A « 2 



-1/2 



dx, 



finds 



Pr(r) 



aT 



-Pi 



{—^j' a = (100 km/s/Mpc x Gyr)" 1 = 9 



78- 



(19) 



(20) 



which is plotted in Figure 3a. Since we neglect the effect of Ao on Ph(h), P t (t) depends on Ao only 
through T = T(fi ,A ). Although T decreases with fi , we numerically found that At = y/ ((r — (t)) 2 } = 
aT\J(h~' 2 ) — (h -1 ) 2 is an increasing function of fio. Nevertheless (r) + At still decreases with fio. In Figure 
3b is plotted the cumulative probability P t> (t) = Ph<{aT / t) that the age of the universe is older than r Gyr. 

Our predictions on the models are summarized in Table 2, which list the expectation values of h and r 
with lcr error, the probability that the Hubble constant is less than 50 km/s/Mpc, and the probability that 
the age of the universe is older than 14 Gyr. 



4. Discussion and Conclusion 



We have presented a method to compute the probability distribution function of the (true) Hubble constant 
and the age of the universe, given the estimate of the Hubble constant in our nearby galaxy samples. Our 
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Figure 2: a) Probability distribution of the Hubble constant Ph(h) (eq.[ll]); b) Cumulative probability Ph<{h) 
that the Hubble constant is less than h (eq.[18]). 




Figure 3: a) Probability distribution of the age of the universe P t (t) (eq.[20]); b) Cumulative probability 
P t> (t) that the age of the universe is older than r Gyr. 
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Table 2: Summary of the prediction on the CDM models 



Model 


{h)±Ah 


(r) ± At 


Ph< 


(0.5) 


Pr> 


(14) 


El 


0.83±0.18 


8.3±1.8 


8.4 x 


: 10" 3 


2.5 x 


10" 3 


E2 


0.83±0.24 


8.6±2.9 


6.1 x 


: lO" 2 


4.1 x 


lO" 2 





0.81±0.10 


10.4±1.4 


3.2 x 


: 10" 4 


1.1 X 


lO" 2 


LI 


0.81±0.10 


13.2±1.7 


3.2 x 


: 10" 4 


3.0 x 


lO" 1 


L2 


0.81±0.18 


13.7±4.0 


4.0 x 


: lO" 2 


3.8 x 


lO" 1 


L3 


0.80±0.17 


13.9±3.9 


4.0 x 


: lO" 2 


3.9 x 


lO" 1 


TCO 


0.80±0.09 


8.3±0.9 


2.2 x 


: 10" 4 


4.2 x 


10" 5 



method takes into account both the observational errors and the cosmic variance, and enables to quantitatively 
compute the constraints on the cosmological models. 

Before summarizing the implications of our results, let us discuss the validity of the spherically symmetric 
approximation which we assumed throughout the present Letter. We described the dynamics of the local 
universe centered around us using the spherical nonlinear model with the mean density contrast smoothed over 
the same effective volume. Since the scales of our interest here are still in linear regime, this approximation 
is statistically valid in practice. To see this in more detail, let us compare our Ph(h) with that computed 
numerically by TCO. Specifically the model TCO in our Figure 2a should be compared with their Figure 
7. Qualitatively they are in agreement, but more careful examination reveals that our PDF is more strongly 
peaked around h = /i j s . In fact we believe that this difference is explained by the fact that they approximated 
the PDF for 8h by the Gaussian; as is clear from their Figure la, the PDF computed directly from the numerical 
simulations is more strongly peaked than the Gaussian fit especially on small scales. In fact our PDF Ph\s{h\6) 
properly takes account of such departure from the Gaussian distribution. Therefore the quantitative difference 
is not due to the assumption of the spherical symmetry. Also it is interesting to note that our model El is quite 
similar to their result based only on data from octant of the sky. TCO ascribed the broadening of the PDF 
to the angular variance in the expansion rate, but Figure 2a implies that most of the effect can be explained 
simply by larger a in small sampling volume (cf, eq.[17]); the smoothing scale R = 12/i _1 Mpc in our model El 
is roughly the same as their effective sampling radius 30/2 /i _1 Mpc. Thus we conclude that our assumption 
of spherical symmetry in the present context is justified on R }t 10 /i _1 Mpc at least in CDM models. 

Applying our method of computing the PDF to the CDM models, one can discuss quantitatively the 
confidence level of the (global) cosmological parameters inferred from the local observations. Let us assume, 
for instance, that the age of the present universe is larger than 14 Gyr. Then as shown in Table 2, the Einstein 
- de Sitter models are rejected with (96 ~ 99.7) % level (models E2 and El) even if we allow for the relatively 
big quoted error in the current local observation. Similarly the model O (Slo = 0.2 and Ao = 0) is viable only 
with 1 % confidence level. Naturally non-zero A models are strongly favored. 

Currently the observational error of the Hubble constant is dominated by the back-to-front geometry of 
the Virgo cluster itself (Freedman et al. 1994), and therefore we do not have strong constraints on the non- 
zero A models. As the number of sample galaxies increases, however, such a systematic error will be reduced 
significantly. Then our probabilistic approach will provide a powerful tool in quantitatively constraining the 
(global) cosmological parameters. 
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